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Abstract 



We present a new determination of fjj a using 5 million e + e — ► cc events 
obtained with the CLEO II detector. Our value is derived from our new 
measured ratio T(D+ -> fj,+v)/T(D+ <j>ir+) = 0.173 ± 0.021 ± 0.031. Using 
B(D+ -> 0vr+) = (3.6 ± 0.9)%, we extract / Ds = (280 ± 17 ± 25 ± 34) MeV. 
We compare this result with various model calculations. 
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I. INTRODUCTION 



Measuring purely leptonic decays of heavy mesons allows the determination of meson 
decay constants, which connect measured quantities, such as the BB mixing ratio, to CKM 
matrix elements. Currently, it is not possible to determine fs experimentally from leptonic 
B decays, so theoretical calculations of fs must be used. Measurements of the Cabibbo- 
favored pseudoscalar decay constants such as fu a provide a check on these calculations and 
help discriminate among different models. 

The decay rate for £>+ is given by |l| @ 

rW-/M = §/W^(i-^)V cs | 2 , (i) 

where Md s is the D s mass, mi is the mass of the final state lepton, V cs is a CKM matrix 
element equal to 0.974 ||, and Gf is the Fermi coupling constant. Various theoretical 
predictions of fo s range from 190 MeV to 350 MeV. Because of helicity suppression, the 
electron mode — > e + v has a very small rate. The relative widths are 10 : 1 : 2 x 10~ 5 
for the t + u, n + v and e + v final states, respectively. Unfortunately the mode with the largest 
branching fraction, t + u, has at least two neutrinos in the final state and is difficult to detect. 

In a previous publication [Q, CLEO reported the measurement of fu a = (344±37±52± 
42) MeV, using the decay sequence D* + — > jDf, — > ji + u. Three other groups have also 
published the observation of — > ji + u and extracted values of fo s - WA75 reported fn s 
as (232 ± 45 ± 20 ± 48) MeV using muons from Df leptonic decays seen in emulsions [|5]; 
BES measured a value of (430i]jQ ± 40) MeV by fully reconstructing mesons close to the 
production threshold in e + e~ collisions ||; and E653 extracted a value of (194 ±35 ±20 ±14) 
MeV from one prong decays into muons seen in an emulsion target 0. 

In this paper we describe an improved CLEO analysis. We use a sample of about 5 
million e + e~ — > cc events collected with the CLEO II detector || at CESR. The integrated 
luminosity is 4.79 fb _1 at the T(4S') resonance or at energies just below. This paper super- 
sedes our previous result which was based on a subset of the current data with 2.13 fb _1 . The 
improvements include a better analysis algorithm, more data, more precise measurements of 
the lepton fakes, and reduced systematic uncertainties. 

II. ANALYSIS METHOD 
A. Overview 

The analysis reported in this paper is based on procedures developed for the previous 
CLEO II measurement of fo s ||. We search for the decay chain D* + — > 7-D+, — > /i + za 
The photon from the D* + decay and the muon from the — > fi + u decay are measured 
directly, while the neutrino is measured indirectly by using the near-hermeticity of the CLEO 
II detector to determine missing momentum and energy. Using the missing momentum as 
the neutrino momentum, we look for a signal in the mass difference 

AM = M( 7 /i+z/) - M(/i + i/), (2) 
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so that the relatively large errors from the missing momentum calculation will mostly cancel. 

To study the AM signal and background shapes and to evaluate the effectiveness of 
our Monte Carlo efficiency simulation, we also collect a data sample of similar topology, 
D*° — > 7-D°, D° — > i^~7r + . We treat these fully reconstructed data events as — > fi + v 
decays by removing the measurements of the tt + from both the tracking chambers and the 
calorimeter to simulate the is, and by "identifying" the K~ as a muon. Our aim here is to 
compare the Monte Carlo simulation of these D*° decays with what we obtain from the data. 

Another useful event sample consists of the decay sequence D* + — > -k + D° , D° — > K~n + , 
since this sample has relatively high statistics and negligible background. We use these events 
to study the missing energy and momentum measurements by eliminating the measurements 
of the fast tc + from the D° decay from both the tracking chambers and calorimeter to simulate 
the neutrino, and call the K~ a muon. 

B. Background 

There are several potential sources of background for this measurement. The real physics 
backgrounds, such as semileptonic decays, are almost identical in muon and electron final 
states because of lepton universality. For the leptonic decay, however, the electronic 
width is negligible in comparison to the muonic width. Thus, performing the identical 
analysis except for selecting electrons rather than muons gives us a quantitative measurement 
of the background level due to real leptons. Df — > ji + u and D + — > ji + is are the only physics 
processes that produce significantly more primary muons than electrons with momenta above 
2 GeV/c in continuum e + e~ annihilations in the T(AS) energy region. D + — > decay 
background in our sample is highly suppressed by the CKM angle (Eq. p, and by the small 
D* + -> jD+ branching ratio, (1.4 ± 0.5 ± 0.6)% §. 

Another source of background results from the misidentification of hadrons as muons 
(fakes). Since muon identification in CLEO II has larger fake rates than electron identi- 
fication, we need to consider the excess fakes in the muon sample relative to the electron 
sample. To determine the hadron-induced muon and electron fake background contributions, 
we multiply the AM distribution of all tracks, excluding identified leptons, by an effective 
hadron-to-lepton fake rate, measured with tagged hadronic track samples. The detailed 
analysis of this effective fake rate is described in Section JTJ. 

After removing the above two components, all remaining events result from either D* + — > 
— ► fits decays, or from spurious combinations of random photons and real — > 
fi + u and D + — ► fi + is decays. The shape of the latter component is determined using the fully 
reconstructed D* + — > n + D°, D° — > K~ir + data sample, and the normalization is determining 
by measuring the D* + / production ratio. Subsequently, we will form a single signal shape 
from these two signal components. 

C. Event Selection and Background Suppression 

Most of the leptons from B meson decays are removed by requiring a minimum lepton 
momentum of 2.4 GeV/c, which is 33% efficient for — > fi + v. Leptons from r + r~ pairs, 
and other QED processes with low multiplicity, are suppressed by requiring that the event 
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either has at least five well reconstructed charged tracks, or at least three charged tracks 
accompanied by at least six neutral energy clusters. To suppress background from particles 
that escape detection at large cos 9, where 9 is the angle with respect to the beam axis, we 
require that the angle between the missing momentum of the event and the beam axis, 6 m i SS , 
does not point along the beam direction, specifically |cos miss \ < 0.9. 

Muons are required to penetrate at least seven interaction lengths of iron, and to have 
|cos 9\ < 0.85. The muon identification efficiency, measured with e + e~ — > fi + [i~ r y events, is 
(85±1)% for muons above 2.4 GeV and is very flat in momentum. Electrons must have an 
energy deposit in the electromagnetic calorimeter close to the fitted track momentum, and a 
dE/dx measurement in the main drift chamber consistent with that expected for electrons. 
The electron identification efficiency for |cos 9\ < 0.85, is found by embedding tracks from 
radiative Bhabha events into hadronic events. For electrons with momentum greater than 
2.4 GeV, a value of (89±2)% is used. 

To subtract the electron data from the muon data we need to have a precise measure of 
the muon to electron normalization. Detector material causes a difference between muons 
and electrons, as electrons tend to radiate more. The correction factor is estimated to lower 
the electron rate by 5%: thus we assign a +5% increase in the electron sample due to this 
outer bremmstrahlung. A Monte Carlo study shows that the main background contributions 
from real leptons in the AM distribution are semileptonic D decays, mostly D — > K£v, niu 
and r\lv. As a specific example of the near equality of the muon and electron rates we made 
a detailed study of the D + — > K £ + u decay. A calculation of the different probabilities that 
a photon is emitted in the decay (inner bremsstrahlung) for D + — > K°l + v was performed 
according to the prescription of Atwood and Marciano [[It]]. This effect raises the electron rate 
by +2.7%. This inner bremsstrahlung correction for the different semileptonic final states 
averages also to +2.7%. We also correct for differences in muon and electron phase space, 
which lowers the relative electron normalization (—1.7% for D + — > K°l + v ). Taking all of 
these sources into account, including the different possible decay modes and the fact that 
the electron detection efficiency is 4% larger than the muon efficiency, we use a correction 
factor of 1.01+0.03 to multiply the electron sample to account for the physics backgrounds 
and the identification efficiency difference. 

Photons must be in the angular region |cos 9\ < 0.71. We require a minimum energy of 
150 MeV, which is 78% efficient for D* + —>■ decay, to eliminate backgrounds caused by 
the large number of low energy photons. Combinations of two photons which have invariant 
masses within two standard deviations of the 7r° mass are eliminated. (The r.m.s. 7r° mass 
resolution is 5 MeV.) We also insist that in the rest frame of the D* + candidate, the cosine 
of the angle between the photon and the D* + direction in the lab be larger than —0.7. A 
small residual b — > ulv background is suppressed by requiring that the thrust axis lines up 
with the D* + candidate momentum so that the cosine of the angle between them is greater 
than 0.975. 



D. Signal Shape and Efficiency 

To evaluate the neutrino four-vector we measure the missing momentum and energy in 
only half of the event; we divide the event into two hemispheres using the thrust axis of the 
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event. The missing momentum p m i SS and energy E miss are calculated using only energy and 
momentum measurements (Ei,pi) in the hemisphere that contains the lepton (kaon). We 
compute the energy sum assuming all tracks are pions, unless they are positively identified 
as kaons, or protons by dE/dx measurement in the drift chamber. We define the missing 
momentum and energy as 

Pmiss = Pthrust ~ J]] Pi and 

■^miss -^beam ^ ] j (3) 



2 

beam 



where the direction of pthrust is given by the thrust axis. The magnitude is p\ hrust = E l 
mj et , where Eb eam is the beam energy and irtj et is the average mass of a charm quark jet, 



measured to be 3.2 GeV using our sample of fully reconstructed D* + events ['I I]. A 
candidate is selected by requiring 1.2 GeV < M(/i + z/) < 3.0 GeV, and that the missing mass 
squared be consistent with a neutrino, | E^ iss — p miss \ < 2 GeV 2 , where the cut values are 
based on studies using the D* + events. Furthermore, we also require p m iss > 0.8 GeV/c 
to suppress backgrounds, since real Df — > fi + i> events must have some missing momentum. 
The D* + candidate momentum is required to be above 2.4 GeV/ c. We find a factor of two 
increase in efficiency by using only one hemisphere to determine the missing momentum 
relative to using the whole event. 

Although the measurement errors on the muon and neutrino tend to cancel when evaluat- 
ing the mass difference in Eq. ^|, the neutrino is poorly enough measured to cause a significant 
broadening of the resolution in comparison with fully reconstructed D* + samples. Improve- 
ment is possible by using the constraint that the muon and neutrino four-vectors must have 
the D s invariant mass. Since the muon is much better measured than the neutrino, we vary 
only the neutrino momentum relative to the selected muon. From conservation of energy 
and momentum, we have 

E Ds =E^ + E U and (4) 



Pd s =Pv+Pv (5) 

Squaring Eq. [| in the local coordinate frame defined by the muon and the reconstructed 
neutrino, using Eq. |] and rearranging shows a relationship between p v and the cosine of the 
angle between the muon and neutrino: 

Pv = (m 2 Ds - mj) / {2E fl - 2p fl cos 9) , where E M = ^/m^ + p 2. (6) 

Fig. |l] shows the constraint as a surface of revolution about the muon momentum vector. 
We start by defining a plane by the vector cross product of the measured muon and neutrino 
three-vectors, though the "correct" solution may lie outside this plane. We next find the 
minimum distance from the measured neutrino momentum vector to the surface. Clearly, 
the new neutrino momentum is the vector sum of the measured neutrino momentum p v meas 
and the distance in momentum space, d, as is shown in Fig. |T|. This procedure improves the 
AM resolution by about 30%. 

We use Monte Carlo simulation to determine the AM signal shape (Eq. |2[) and to estimate 
our efficiency. Since this analysis involves reconstructing a missing neutrino, we are concerned 
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that the Monte Carlo will not adequately simulate the data. As a check we evaluate the 
accuracy of our simulation using our D*° — > 7-D , D° — > K~n + sample, where we eliminate 
the 7r + to simulate the neutrino and treat the K~ as a muon. 

We start with a D*° — > 7-D , -D° — > i^~7r + Monte Carlo simulation. Fig. 0(a) shows the 
fully reconstructed mass difference AM = M(^Kir) —M(Kir) distribution after a cut on the 
K~ir + invariant mass of ±30 MeV around the known D° mass (where the r.m.s resolution 
is 8 MeV). The kaon is required to have momentum greater than 2.4 GeV/c, which is the 
same cut as we use on the muon in the —>■ channel. In the AM distribution there 
is a substantial signal but also significant background, so a AM sideband subtraction is 
performed. We use a bin-by-bin subtraction. The central value of the signal is 142 MeV 
and the r.m.s. width is 5.5 MeV. The sidebands used are 114-126 MeV and 159-170 MeV. 
After applying the additional background suppression cuts, described above, we obtain the 
mass difference distribution AM = M{^Kp m i SS ) — M(Kp m i SS ) shown in Fig. |3|. There is a 
clear signal peak associated with the photon and it is fitted to an asymmetric Gaussian with 
low side and high side it's of 15 MeV and 16 MeV, respectively. The small flat component 
results from replacing the correct photon with another photon. 

The partial efficiency for neutrino detection only from Monte Carlo for a fully recon- 
structed D*° event with both the D*° and its kaon daughter having a momentum greater 
than 2.4 GeV to appear in the signal peak after neutrino reconstruction is found to be e„ 



= (38.9±2.6)% 0. The overall detection efficiency for D*° -> D° 7 , D° -> K + p miss is 
(4.8 ± 0.3)% 0. 

Next, we repeated the analysis described above for the fully reconstructed D*° — > 7.D 
data sample. The fully reconstructed AM distribution is shown in Fig. 0(b). The AM 
distribution for the missing neutrino is shown in Fig. f| where the sideband subtraction again 
has been performed. The fitting parameters derived from the Monte Carlo signal shape fit 
the data very well, with a x 2 of 23 for 27 degrees of freedom and a confidence level of 69%. 
The partial efficiency for neutrino detection oly of (38.5±3.7)% agrees well with Monte Carlo 
simulation. 

In principle the resolution and efficiency for D* + — ► 7D+, D+ — > fi + u can be somewhat 
different from that for the D*° sample described above, because of the different fragmentation 
with an s quark rather than a u quark. Since our Monte Carlo simulation accurately describes 
the D*° — > jD°, D° — > K~it + process, we rely on it for our D* study. In Fig. we show the 
AM = M(7/x + z/) — M(/i + z/) distribution from Monte Carlo simulation. This distribution 
contains a Gaussian part due to the signal, plus a background which occurs when the correct 
photon from the D* decay is replaced with another random photon in the event. We fit 
the histogram with an asymmetric Gaussian signal shape having low side and high side cr's 
of 15 MeV and 17 MeV, and the function y/x — xoe~ a<yX ~ x ^ to parameterize the random 
photon component, where x = AM. The Gaussian signal shape agrees well with the D*° 
Monte Carlo and data. Using the Gaussian signal part only, the overall efficiency is found 
to be (4.2±0.3)%, where the error includes the systematic effect of the efficiency difference 
between data and Monte Carlo determined by the D*° sample. 

An additional source of background in the AM distribution comes from direct — > 
decays which pair with a random photon to form a D* + candidate. These are in addition 
to D* + events where the correct photon is replaced by another photon, as mentioned above. 
These two contributions are fixed relative to the direct D* + — * jD+, \w signal using 



our measurement of D* + /D+ production ratio above 2.4 GeV of 1.08±0.13 (see below). 
Thirdly, there is a small contribution from D + — > fi + u decays combined with a random 
photon. The shape in AM of all these contributions is modeled using the D* + — > tt + D° 
event sample, by combining the M(Kp m i SS ) candidates with random photons in the same 
event, and fitting with the functional form yjx — XQe~ a( - x ~ x °^ to parameterize the total random 
photon component. The distributions in Fig. [| and the random photon component function 
are summed using appropriate weights to produce the expected shape for the sum of the 
D* + — ► 7.D+, — > fi + i> signal plus random photon background shown in Fig. |6|. 



E. Measurement of the and D* + Rates 

In order to measure the relative rates of D* and D s production, and the absolute level of 
D s production above 2.4 GeV/c we use the — > </>7r + decay mode. The <p is searched for 
in the K + K~ decay mode. We require the photon from the D* + decay to satisfy the same 
requirements as for the fi + i> final state. The detection efficiency for the decay mode is 



22.3%, while for the D*+ the efficiency is 9.4% |4 . 

Fig. |7|(a) shows both the invariant mass of the 07r + . In (b), we show AM = M(707r + ) — 
M(07r + ) after requiring that the (j)ir + mass be within ±24 MeV of the mass. 

Fitting the data to Gaussian signals shapes whose widths are determined by Monte 
Carlo simulation we find 5728±123 — ► 07r + events and 1256±54 D* + events. Taking 
into account the relative efficiencies we determine that the ratio of D* + /D+ production is 
1.08±0.13. This number reflects the direct production of a vector charmed- strange meson 
relative to the direct production of a pseudoscalar charmed-strange meson, above 2.4 GeV/ c 

Ga- 



in. LEPTON FAKE BACKGROUND CALCULATION 

Even after strict lepton identification requirements have been applied, significant numbers 
of hadron fakes still enter our signal region because of the abundance of fast hadron tracks. To 
properly account for the hadron fake background, we need to measure precisely the effective 
excess muon to electron fake rate ratio to derive the correct background level. The D* 
decays provide us with well-tagged kaon and pion samples. In our previous publication, the 
uncertainty in the fake rate value dominated the systematic errors. One major improvement 
of the current analysis is the better determination of these rates for muons and electrons from 
much larger tagged data samples obtained by using new data and adding more channels. 

In this analysis, in addition to the decay sequence D* + — > D°tt + — > (K~tt + )tt + , we also 
include D* + -> D°n + -> (K-p + )-K + , and D*° -> D n -> (K-tt + )tt to get as many events 
as possible. Kg — ► 7r + 7r~ samples are also used to determine the pion fake rate and are 
combined with the D* results to get better statistics. Over 10,000 events were collected with 
either a tt or K with momentum greater than 2.4 GeV from the above channels. 

In Fig. we show the M{ir + K~-K + ) — M(K~tt + ) mass difference after a cut on Kir mass 
consistent with the D° mass for kaons or pions which pass our cuts for muons or electrons 
The number of events is determined by a fit with a double Gaussian for the signal and half- 
integer power polynomials for background. Both fitting function shapes are derived from 
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the mass difference distribution without lepton identification suppression. Our extracted 
fake rates (before decay in flight correction) are listed in Table |. The same reconstruction 
methods are used to collect kaon and pion samples from the channels D° — > K~p + , and 
D*° — > D 7T° where D° goes to K~tt + sample. The fake rates are determined by fitting 
the mass distributions for the amount of signal. The fake rates derived from the different 
channels are summarized in Table |, and the weighted average fake rates are also shown. 



TABLE I. Fake Rates for P > 2.4 GeV/c 



Data Samples 


#of 






Fake Rates (%) 






K 


IT 


K/fi 


K/e 


7T/V 


n/e 


D*+{D° -» K--K+) 


9404 


7461 


0.94±0.11 


0.04±0.05 


0.60±0.12 


0.24±0.06 


D* + (D° -» K~p+) 


1368 


682 


1.23±0.33 


0.22±0.20 


0.30±0.40 


0.15±0.21 


D*°(D° K--K+) 


3174 


2048 


1.07±0.21 


0.17±0.10 


0.84±0.35 


0.60±0.31 


Ks — ► 7T + 7T~ 




3527 






0.74±0.15 


0.37±0.10 


Total /Aver age 


13964 


13718 


0.98±0.08 


0.12±0.05 


0.65±0.08 


0.31±0.06 



The contributions to the lepton fake rates from kaon and pion decays in flight are not 
necessarily included in the above procedure because particles decaying close to the production 
point may not appear in the D° mass peak. To account for this effect, we used a Monte 
Carlo study of 200,000 D* + — > D°tt + — > (K~tt + )tt + events. After muon identification cuts 
are applied, the D* + mass difference plot has a peak region used to derive the fake signal 
and a tail away from the peak, which is due to events in which the kaon decays. We extract 
a correction factor to the fake rate of 1.18±0.06 by computing the ratio of the tail area to 
the peak area. We find no events out of the D* mass difference peak in which the pion 
has decayed. This is because of the relatively long pion lifetime and because the muon 
momentum is very close to that of the parent pion. 

We determine the hadron induced muon and electron fake background contributions by 
multiplying the AM distribution of all tracks, excluding leptons, by the effective fakes rates 
determined above. The fractions of kaon, pions and protons are 67%, 20% and 13% as 
ascertained from Monte Carlo simulation. The effective Fake rates from protons and anti- 
protons are small, ~0.1%, and almost equal for muons and electrons. 

IV. RESULTS 

The AM distributions for the muon and electron data and the calculated effective excess 
of muon fakes over electron fakes are given in Fig. f|(a). The histogram is the result of 
a x 2 fit °f the muon spectrum to the sum of three contributions: the signal distribution 
evaluated with the D* + Monte Carlo plus random photon background evaluated with the 
D* + sample, the scaled electrons, and the excess of muon over electron fakes. Here, the sizes 
of the electron and fake contributions are fixed and only the signal normalization is allowed 
to vary. 

We remind the reader that the signal consists of two components, whose relative nor- 
malization is fixed. These two components are the decay D* + — > jD£, —> fi + u and the 
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direct decays Df — > fiu and D + — > fi + u. Our measurement of the D* + /D+ production ratio 
allows us to constrain the relative normalization. We find a signal of 182±22 events in the 
peak which are attributed to the process D* + — > jDf, — > ji + u. We also find 250±38 
events in the flat part of the distribution corresponding to — > ji + u or D + — >• jU + z/ decays 
coupled with a random photon. The contribution of a real D + — > fi + u decay with random 
photons is not entirely negligible since the D* + — > branching ratio does not enter. The 
D + fraction is estimated to be about (18±8)% relative to the total — > [i + v plus random 
photon contribution. 

To explicitly display the signal, we show, in Fig. |](b), the AM distribution after the 
electrons and the fakes are subtracted. The curve is a fit of the data in Fig. £|(a) to the 
signal shape calculated from the D* + sample and random photon background calculated 
from the D* + sample. All of the events in this plot are signal, the background having 
already been subtracted. 

Using the fit result of 182±22 events, we extract a width for £)+ — > by normalizing to 
the the efficiency corrected number of fully reconstructed D* + — > 7-D+, — > 07r + events, 
HI 



24740±1200±810 



Monte Carlo. We find 



The efficiency for reconstructing the (f)ir + decay is obtained from 



0.173 ±0.023 ±0.035 



)7T 



(7) 



where the first error is the statistical error on the measured numbers of /z + z/ and 07r + events. 
The second error is the total systematic error of 18%, whose components are summarized in 
Table 0. 

TABLE II. Systematic Errors on Width Ratio (%) 



Source of Error 



Value 



Size of error 



Muon fake rate 

Electron fake rate 

ir/K/p fractions (sources of fakes) 

fi/e normalization 

Detection efficiency 

D* + /Df production ratio 

d)TT + normalization 



(0.69±0.05)% 
(0.21±0.03)% 
67%/20%/13% 
1.01±0.03 
(4.2±0.3)% 
1.08±0.13 
24740±1200±810 



9 
7 
7 
9 
7 



Total systematic error 



20 



The errors that arise from the relative muon to electron normalization, the muon fake 
rate, the electron fake rate, and the D* + /D+ production ratio, are estimated by fitting the 
data with each parameter changed by ± la. The error on the relative fractions of pions, kaons 
and protons entering into the fake rate calculation is computed by changing the fractions to 
70%, 20% and 10%, respectively. We judge this to be the outer limit at 90% confidence level 
of the change possible in these ratios. This, in turn, changes the excess muon to electron to 
fake rate by 12% leading to a 7% change in the yield. A systematic error of ±3% for the 
detection efficiency of the normalization mode 07r + is also included. 
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The radiative decay rates for Df — ► £ + i"j and B + — > £ + z^7 have been considered by 
Burdman, Goldman and Wyler ||23|| . They predict that 



r(D+ - 

where vector coupling constant which has a value of approximately 0.1 for the Df 

meson. While the radiative decay rate for B + is comparable to the non-radiative rate, the 
radiative decay rate for Df is estimated to be between 0.1% and 1% of the non-radiative 
rate. Furthermore, they also predict that the radiative muon and electron rates are equal, 
so our electron subtraction would remove any residual effect. 



V. CONCLUSIONS 

We have measured the ratio of decay widths 

F(D+ -f /i + z/)/r(D+ -f 07T+) = 0.173 ±0.023 ±0.035 . (9) 

To extract the decay constant fn s we need to known the partial width for the — > <f)iT + 
decay. The total width is well known because of precise lifetime measurements ||, but 
the absolute 07r + branching ratio has a large error. Using the latest PDG average value 
B(D+ -> 07T+) of (3.6±0.9)%, and t Ds = (4.67 ± 0.17) x 10~ 13 s, we find 

f Ds = (280 ± 19 ± 28 ± 34) MeV . (10) 

The first error is statistical, and the second is systematic resulting from our relative width 
ratio measurement, and the third error reflects the uncertainty in the absolute Df — > 07r + 
branching ratio. This result supersedes our previous one, using a data sample that includes 
the one used in the previous analysis. The reduction in the central value is primarily due to 
the better measurement of the lepton fake rates that lowered the pion/electron fake rate. 



For comparison, we list in Table |1| the old CLEO result and published results from 
other experiments that used the D+ — > /i + z/ decay to measure fc 3 . We have changed the 
values of fr> a according to the new PDG D s decay branching fractions for the normalization 
modes, and have corrected the old CLEO result by using the new fake rates determined in 
this analysis ||. The lowering of the central value of the old CLEO result is mostly due to 
the change in the fake rate determination, which is now much more precise. 

In addition, there are new results using the — > r + v decay from the L3 collaboration 
IH of 309 ± 58 ± 33 ± 38 MeV, and 330 ± 95 from the DELPHI collaboration g?|. Our new 
measurement gives the most accurate of fo a - 

Theoretical predictions of fjj a have been made using many methods. Recent lattice 



gauge calculations pB| give central values of 199 to 221 MeV with quoted errors in the ±40 



MeV range. Other theoretical estimates use potential models whose values |L9[ range from 
210 to 356 MeV, and QCD sum rule estimates [0 that are between 200 and 290 MeV. 
Predictions for fu a have also been made by combining theory with experimental input. 
Assuming factorization for B — > D*D~ decays combined with measured branching ratios, 



gives a value of fo a range of about 280 MeV with an error of about 60 MeV pTJ. Use of 
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TABLE III. Current Experimental Results on fjj g Using Df — > [i + v 



Collaboration 


Observed 


Published fjj a 


Corrected fu g 




Events 


value (MeV) 


value (MeV) 


CLEO (old) § 


39±8 


344 ± 37 ± 52 ± 42 


282 ± 30 ± 43 ± 34 


WA75 § 


6 


232 ± 45 ± 20 ± 48 


238 ± 47 ± 21 ± 48 


BES § 


3 


430l^g ± 40 


Same 


E653 § 


23.2 ± 6.0±g;g 


194 ± 35 ± 20 ± 14 


190 ± 34 ± 20 ± 26 


CLEO (this work) 


182±22 




280 ± 19 ± 28 ± 34 



experimental data on isospin mass splittings in the D* and D system gives a value for fu of 
290 MeV pH. (f Ds is thought to be 10% to 20% higher than f D ). 
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FIG. 2. The AM = M(7i-C~7r + ) - M(K~ir + ) mass difference distributions for fully recon- 
structed events of D*° — > 7D , D° — ► i"C~7r + , after a requirement that the K~tt + mass be within 
2.5 standard deviations of the D° mass, (a) -D*° Monte Carlo (b) data. 
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FIG. 3. The AM = M(~fKp m i ss ) — M(Kp m i ss ) mass difference distributions for the simulated 
missing momentum analysis using Monte Carlo of D*° — ► 7D , D° — > i^ _ 7r + . A sideband sub- 
traction to remove background in the initial D*° selection has been applied. The curve and fitting 
procedure are described in the text. 
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FIG. 4. The AM = M(^Kp m i ss ) — M(Kp m i ss ) mass difference distributions for the simulated 
missing momentum analysis for the D*° — > 7D , D° — > K~ir + real data. A sideband subtraction to 
remove background in the initial D° selection has been applied. The curve used is from the Monte 
Carlo signal shape. 
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FIG. 5. The AM = M (7 UPmiss) — M(np m i ss ) mass difference distributions for the missing 
momentum analysis for the D* + — > jDf, Df — > ;U + ^ Monte Carlo. The curve and fitting procedure 
are described in the text. 
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FIG. 6. The D* + — > 7-D+, Df fj, + u signal distribution plus random photon background as 
determined from the signal Monte Carlo simulation combined with the D* + data sample analyzed 
for the missing v as D*° — > 7-D , D° — > K~p m i ss . The curve is a fit using the functions described 
in the text. 
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consistent with the known Df mass. The signal 
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FIG. 8. M(irKir) — M(Kir) mass difference distributions for four cases of hadrons identified 
as leptons: (a) kaon as muon, (b) kaon as electron, (c) pion as muon, (d) pion as electron. The 
signal shapes were determined from the distribution of mass difference for fully reconstructed D* + 
candiates. The area of the Gaussian component and the normalization of the background are 
allowed to float. 
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FIG. 9. (a) The AM mass difference distribution for D* + candidates for both the muon data 
(solid points), the electron data (dashed histogram) and the excess of muon fakes over electron 
fakes (shaded). The histogram is the result of the fit described in the text, (b) The AM mass 
difference distribution for candidates with electrons and excess muon fakes subtracted. The 
curve is a fit to the signal shape described in the text. 
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